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Abstract 


In order to be cost-effective, space structures must be extremely 
light-weight, and subsequently, very flexible structures. The power 
system for Space Station Freedom is such a structure. Each array 
consists of a deployable truss mast and a split "blanket*' of photo- 
voltaic solar collectors. The solar arrays are deployed in orbit, and 
the blanket is stretched into position as the mast is extended. 
Geometric stiffness due to the preload make this an interesting non- 
linear problem. 

The space station will be subjected to various dynamic loads, 
during shuttle docking, solar tracking, attitude adjustment, etc. 
Accurate prediction of the natural frequencies and mode shapes of the 
space station components, including the solar arrays, is critical for 
determining the structural adequacy of the components, and for 
designing a dynamic controls system. 

m^-paper chronicles the process used in developing and 
verifying the finite element dynamic model of the photo-voltaic 
arrays / Various problems were i denti f ied„in the^ investigation , such 
as grounding effects due to geometric stiffness, large displacement 
e f fgjjt s, and pseudo-stiffness (grounding) due to lack of required 
rigid! body modes. Var i ou s analysis techniques, such as development of 
rigorous solutions using continuum mechanics, finite element solution 
sequence altering, equivalent systems using a curvature basis, Craig- 
Bampton superelement approach, and modal ordering schemes were 
utilized. This — paper emphasizes the grounding problems associated 
with the geometric stiffness :> & 


Nomenclature 


factor defined by Eq.(13) 
arbitrary constants in Eq. (10) 


1 


d/dx, or 
d/dt, or 
E 

e a 

(F) 

F(x.t) 

g 

I 

[K] 

[«e] 

[ K el 

L 

M 

dM 

m 

P 

P ' 

(R) 

T 

<u) 

u 

UA 

U B 

v 

V 

dV 

V 

dVol 

W 

X 

y 

P 

6 

6 

e 

o 


differential operator with respect to position 
differential operator with respect to time 
modulus of elasticity 
axial strain 

input force vector at the beginning of a step 

applied transverse force 

factor defined by Eq.(14) 

moment of inertia 

stiffness matrix 

elastic stiffness matrix 

geometric stiffness matrix 

length 

moment 

change in moment 
mass per unit length 
axial force 

pseudo-force necessary for equilibrium 

force vector, output force vector at the end of a step 
kinetic energy 

displacements at the node points 
longitudinal displacement 
strain energy due to axial load 
strain energy due to bending 
transverse displacement 
shear 

change in shear 

potential of the external loads 

change in volume 

natural frequency 

axis defined by Figure 1 

axis defined by Figure 1 

1/2 the angle of rotation 

factor defined by Eq.(ll) 

factor defined by Eq.(12) 

angle of rotation 

stress 


Introduction 


NASA's Space Station Freedom consists of various modules 
supported by a space truss. Power for the space station will be 
provided by a deployable system of split blanket phot o— vo 1 t a i c arrays, 
which will have two degree of freedom rotational capabilities in order 
to track the sun during its orbit. The arrays are designed to be 
operated in a zero-gravity environment. 

NASA Lewis Research Center, along with its contractors, have the 
responsibility for developing a verified finite element dynamics model 
of the solar arrays, which could be combined with the other space 
station substructures for both structural and dynamic control studies. 
The development of the model necessitated the use of unique 
procedures, and rigorous analytical checks. 

The procedure included the following: 
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1. Development of an idealized model of the solar arrays, and 
derivation of a unique solution for the response frequencies for 
the idealized array cantilevered from the space truss, using 
equations developed from continuum mechanics .[ 1 ] 

2. Comparison of the frequencies from the MSC/NASTRAN finite element 
dynamic model of the idealized array with the rigorous solution 
from continuum mechani cs . [ 2 ] 

3. Refinement of the finite element mesh. 

4. Rigid body mode checks of the finite element models. 

5. Various parameter studies involving the amount of tension in the 
blanket, rigidity of the blanket tip beam, type of elements used, 
etc . . 

6. Cra i g-Bampton approach for appending rigid body modes to 
substructures (superelements) [3], 

7. Modal ordering schemes for identifying "important" modes. 

8. Study of grounding effects due to lack of rigid body mode 
capabilities. [4] 

A detailed summary of the project was presented [5]. It should be 

noted that this study is ongoing at the present time. 

This paper will be restricted to the grounding problems associated 

with the geometric stiffness due to blanket pre-load. 


Grounding 

The space station solar arrays were modeled utilizing 
MSC/NASTRAN. As a routine check, the stiffness matrices generated by 
the model were multiplied by a matrix of rigid body modes, and large 
p s e u d o - f o r c e s were developed (grounding) . The cause of this 

"grounding" phenomenum was examined. 

Finite element solves non-linear problems of the form 


[[Ke] + [Kg]] * {u} - { R } - {F> 


where [Ke] is the elastic stiffness matrix, and [Kg] is the geometric, 
or initial stress stiffness matrix. 

[Kg] is a function of the pre-load. Thus, it equals zero for a 
linear problem. [Ke] possesses the required rigid body modes. 
However, [Kg] lacks the capacity for rigid body rotation. Hence, an 
erroneous stiffening, or "grounding", occurs when a pre-loaded beam 
deforms . 


3 



matrix , 


The traditional, or consistent geometric stiffness 
developed by Martin [6] and others, is 
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This matrix does not possess rigid body rotation capabilities. 
Various refinements to the geometric stiffness have been developed 
which contain higher order terms [6,7,8] . However, none of these 
possess all the rigid body modes. Bosela [4] developed a modified 
[Kg] with complete rigid body modes when used with an exact rigid body 
rotation matrix, but [Kg] lost some of its rigid body capabilities. 


Closer examination of the traditional formulation of [Kg] 
indicated that there is a load imbalance in the representation, and 
that pseudo-forces occur to maintain equilibrium. (Fig. 1) 



Fig 1. P’ Represents Pseudo- f or ces Required for Equilibrium 

In Reference [9], Collar and Simpson indicate that the lack of 
rigid body rotation capabilities for [Kg] is not a problem, because 
the energy representation is correct. It can be shown that 4 it is 
correct to p 2 terms, but error does occur, as a function of p . For 
large rigid body rotation, as will occur with the solar arrays, this 
is significant. 

It should be noted that as long as the pre-load P is assumed to 
remain horizontal during rotation, work will be done by the force. 
Thus, true rigid body rotation cannot occur. In order for the strain 
energy to equal zero, the force P must change its orientation as the 
beam rotates ( ie. a follower force). 


Rigorous Solution Of Pre-Loaded Beam 

Suppose we have an axially loaded beam in space subjected to a 
time varying transverse loading (Figure 2). 
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Figure 2 Beam in Tension and Differential Element 


The kinetic energy is 
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The strain energy due to bending is 
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The strain energy due to axial load is 


U 


A 


1 

2 


o e„ dVol 

a 


Letting dVol * dA dx and applying non-linear elasticity yields 
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Neglecting axial displacement and higher order terms yields 


( 4 ) 
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The potential of the external loads is 


V 


F ( x , t ) v dx + V 0 v ( 0 , t ) + M 0 v'(O.t) 
- V L v ( L , t ) - M l v’ (L.t) 


(6) 


Applying Hamilton's principle, and performing the variation, yields 


1 2 L 



EIv"6(v")+Pv'6(v')-mv6(v)-F(x,t)6(v)Jdx 

+V o 6v(0,t)+M o 6v' (0,t)-V L 6v(L,t)-M L 6v' (L.t)jdt 

Integrating by parts yields the differential equation 


tj 0 


d* /dx* ( El d * v/dx * ) - P d*v/dx* + m d*v/dt* = F(x,t) 


(7) 


(8) 


which agrees with Clough in reference [10], after a sign change 
required to express the axial force in tension instead of compression. 
This is also in agreement with Shaker in Reference [11]. 

For a beam in space, the moment and shear at the end points must 
equal zero. Thus, the boundary conditions are 


EIv”(0,t)=EIv"(L,t)=v'"(0,t)-P v'(0,t)=v"'(L,t)-P v’(L,t)=0 

El El 


(9) 


Choose a solution of the form 


v(x) = Djsinffix) + D 2 cos(6x) + D 3 sinh(Gx) + D 4 cosh(ex) 

1/2 

where 6 = |( a 4 + g 4 /4 ) ^ ^-g* / 2 
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g 2 = P/EI 


(14) 


Applying the boundary conditions at x=0, and after much mathematical 
manipulation, yields 


v(x) = D 3 

6sin6x + sinhEx 

e 

+ °4 

€ 2 cos6x + coshEx 
6 2 


m m 
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(15) 


Applying the boundary conditions at x-L, and after more mathematical 
manipulations, yields 


r " 

D 3 6 3 cosh6L - 6 3 cos6L 


€ 3 s i n6L + 6 3 sinhEL 


(16) 


Express ing E q . ( 1 5 ) and Eq.(16) into matrix form, setting the 
determinant equal to zero, and after more mathematical manipulations, 
the following characteristic equation is obtained 


±2a 6 ( coshGLcos6L-l ) + ( £ 6 - 6 6 ) s i nhEL s i n6L = 0 


(17) 


Using Eq . (13) , this can be expressed as 

±w 3 (m/EI ) 3/2 (coshGLcos6L-l ) + ( £ 6 -6 6 ) s i nhELs in6L = 0 . (18) 

By observation, when w=0, a=0, and 6=0. Letting sin(0)=0 yields 

w 3 (m/EI) 3/2 (coshELcos6L-l) = 0 . (19) 

The w 3 term indicates that there must be three zero roots of 
"w", which suggests the three required rigid body modes. 

Conclusion 


Lack of complete rigid body mode capabilities is inherent in the 
physical representation of the pre-tensioned beam problem currently 
used to formulate the geometric stiffness matrix. This lack of 
complete rigid body mode capabilities invalidates the rigid body mode 
check for non-linear problems, and adversely impacts the use of 
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traditional finite element techniques to predict dynamic response of 
pre-loaded structures unless the missing rigid body modes are somehow 
apppended on to the structure, such as by the Craig-Bampton technique. 

The rigorous solution of the ax i al ly- 1 oaded beam with free/free 
boundary conditions developed in this paper may lend itself to the 
development of a new geometric stiffness matrix for a beam element 
with full rigid body capabilities. 
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